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We discuss a consistent treatment of the light-front gauge-boson and meson wave functions in 
the analyses of the generalized parton distributions(GPDs) and the scattering amplitudes in deeply 
virtual Compton scattering(DVCS) for the pion. The continuity of the GPDs at the crossover, where 
the longitudinal momentum fraction of the probed quark is same with the skewedness parameter, 
and the finiteness of the DVCS amplitude are ensured if the same light-front radial wave function as 

■ that of the meson bound state wave function is used for the gauge boson bound state arising from 
the pair-creation(or nonvalence) diagram. The frame-independence of our model calculation is also 

f"^ , guaranteed by the constraint from the sum rule between the GPDs and the form factors. 

(N ' PACS number(s): 13.40.Gp, 12.39.Ki, 13.60.Fz 

3 : 

I. INTRODUCTION 

o\ : 

Compton scattering provides a unique tool for studying hadronic structure. The Compton amplitude probes the 

■ hadrons through a coupling of two electromagnetic currents and thus can be regarded as a generalization of hadron form 
' factors. For the Compton amplitude, when the initial photon is highly virtual but the final one is real, one arrives at 

the kinematics of deeply virtual Compton scattering(DVCS). According to the Quantum ChromoDynamics(QCD), the 
amplitudes at large momentum transfer factorize in the form of a convolution of a hard scattering amplitude which can 
be computed perturbatively from quark-gluon subprocesses multiplied by process- independent distribution amplitudes 
containing the bound-state nonperturbative dynamics for each of the interacting hadrons. Thus, the most important 
contribution to DVCS amplitude is given by the convolution of a hard quark propagator and a nonperturbative function 
describing long-distance dynamics which is known as "generalized parton distributions(GPDs)" and also known as 
"skewed parton distributions(SPDs)" The latter has been studied extensively in searching for possible new 

physics(see, for example, Ref. [Q and references therein). 

The GPDs serve as a generalization of the ordinary (forward) parton distributions and provide much more direct 
Oh, and sensitive information on the light-front (LF) wave function of a target hadron than the hadron form factors. In 
particular, the momentum of the "probed quark" in GPDs is not integrated over, but rather kept fixed at longitudinal 
momentum fraction x, while for the form factor it is integrated out (due to the nonlocal current operator ■(/'(O^^i/K 2 ) 
for the GPDs in contrast to the local vertex ?/5(0)7 M ?/>(0) for the form factor). The form factors are then just moments 
of the GPDs. At the expense of being generalized amplitudes, the GPDs always involve the nonvalence contributions 
due to the nature of longitudinal asymmetry £ = (P — P') + /P + (so called "skewedness" parameter) between initial(P) 
and final(P') hadron state momenta. While the kinematic region where the longitudinal momentum fraction x of the 
probed quark is greater than the skewedness parameter £ (i.e. 1 > x > £) is called "DGLAP region" ||, the rest of the 
longitudinal momentum region < x < ( is called "ERBL region" ||] . The DGLAP and ERBL regions have also been 
denoted as the valence and nonvalence regions in the LF dynamics, respectively, because the parton-number-changing 
nonvalence Fock-state contributions cannot be avoided for < x < C while only the parton-number-conserving valence 
Fock-state contributions are needed for 1 > x > (. Thus, it has been a great challenge to calculate the nonvalence 
contributions to the GPDs in the framework of LF quantization. 

Although many recent theoretical endeavors [f7|-[l2"[ have been made in describing the GPDs in terms of LF wave 
functions, the task has not yet been satisfactory enough for practical calculations. In Refs. |8| and ||, the nonvalence 
contributions to the GPDs have been rewritten in terms of LF wave functions with different parton configurations. 
However, the representation given in Refs. Q and [|| requires to find all the higher Fock-state wave functions while 
there has been relatively little progress in computing the basic wave functions of hadrons from first principles. In 
Refs. []l(| and jll|, the GPDs were expressed in terms of LF wave function but only within toy models such as the 
't Hooft model of (l + l)-dimensional QCD [jl0| and the scalar Wick-Cutkosky model |ll|, respectively. While these toy 
model analyses are helpful to gain some physical insight on the properties of the GPDs (especially, the time reversal 
invariance, the continuity at the crossover between the DGLAP and ERBL regions, and the sum rule constrained by 
the electromagnetic form factor), the real (3 + l)-dimensional QCD motivates us to come up with the more realistic 
model for the application to the analysis of GPDs. 
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In an effort toward this direction, we have presented an effective treatment of handling the nonvalence contributions 
to the GPDs of the pion |Q using our LF constituent quark model(LFQM), which has been phenomenologically quite 
successful in describing the spacelike form factors for the electromagnetic and radiative decays of pseudoscalar and 
vector mesons ^| ^) and the timelike weak form factors for exclusive semileptonic and rare decays of pseudoscalar 
mesons ]l6[-^8[ . Our effective treatment of handling the nonvalence contributions is based on the covariant Bethe- 
Salpeter(BS) approach formulated in the LF quantization flrj|| which we call LFBS approach and has been previously 
applied to the exclusive semileptonic and rare decays of pseudoscalar mesons 0Jl^] providing reasonable results 
compared to the data. 

However, an artifact of discontinuity at x = C occured in our previous calculation of GPDs [Q. As we provided 
the reasoning fl2]| , the discontinuity is caused by the different x = £ behavior between the gauge boson vertex and 
the hadronic vertex if the wave function for the gauge boson vertex is taken differently from that for the hadronic 
vertex. The similar observation was made recently in Ref. jflj using a different model(scalar Wick-Cutkosky model). 
Such discontinuity at x — Q may cause a divergence in the DVCS amplitude. 

In this work, we improve our previous analysis p2| by taking the same LFBS approach for both vertices of meson 
and gauge boson and explicitly show that our reasoning presented in Ref. is correct; i.e. the continuity of GPDs 
at the crossover(a; = Q is ensured by this consistent treatment of vertices. We also discuss the value of GPDs at 
x = £ in conjunction with the single spin asymmetry (SSA) fl9| , ^o|| and calculate the scattering amplitude contributing 
dominantly to the Compton scattering of the pion in the deeply virtual region. The paper is organized as follows. In 
Section II, we derive the GPDs as a nonperturbative formulation of the light-front dominated deeply virtual Compton 
scattering(DVCS) 7*7r — > fir of the pion and introduce the necessary kinematics following the notation employed by 
Radyushkin In Section III, we discuss the LFBS approach and present the consistent treatment of the gauge 
boson and hadron vertex functions in handling the nonvalence contributions to the GPDs. The implication of the 
vector meson dominance(VMD) at the gauge boson vertex in this improved analysis and the GPD value at x = C are 
also discussed. In Section IV, we show our numerical results for the GPDs of the pion that satisfy the continuity at 
the crossover and in turn give the finite scattering amplitude for DVCS of the pion. The frame-independence of our 
model is checked by the sum rule between the GPDs and the pion form factor. We also comment on the polynomiality 
conditions associated with the D-term contribution plLp| in ERBL region. Conclusions follow in Section V. 



II. SCATTERING AMPLITUDE IN DVCS 




P P' P P' 

1 (s-channel) 1 1 (u— channel) 

FIG. 1. Handbag diagrams contributing dominantly to Compton scattering in the deeply virtual region. The lower soft part 

consists of a hadronic matrix element which is parametrized in the form of generalized parton distribution functions. 
We begin with the kinematics of the virtual Compton scattering (see Fig. [l]) of the pion 

7 *( (? )+ 7r (P)^ 7 ( ( /)+7r(P'), (1) 

where the initial (final) hadron state is characterized by the momentum P (P') and the incoming spacelike virtual and 
outgoing real photon momenta by q and q', respectively^ We shall use the component notation V = (V + , V~ , Vj_) 
and our metric is specified by V r± = (V° ± V 3 ) and V ■ V = V + V~ — Vj_. 



1 If both photons are far off-shell and have equal spacelike virtuality, it becomes virtual forward Compton amplitude and 
its imaginary part determines structure functions of deep inelastic scattering(DIS). If they are both real but the momentum 
transfer is large, it becomes wide-angle Compton scattering(WACS) amplitude 0|. 
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Defining the four momentum transfer A = P-P', one has P = [P+, M 2 /P+, Q±], P' = [{1-QP+, (M 2 + A\)/(l- 
C)P + ,-A_l], and A = P - P' = [(P+,{A 2 + A 2 L )/CP + , A±], where M is the pion mass and C = A+/P+ is the 
skewedness parameter describing the asymmetry in plus momentum. The squared momentum transfer then reads 



* = A*=2P.A=-^±^i. (2) 

i -c 

Since A^ > 0, t has a minimum value — t min = Q 2 M 2 /(l — £) at given £. As shown in Fig. [I], the parton emitted by 
the pion has the momentum k, and the one absorbed has the momentum k' = k — A. As in the case of spacelike form 
factors, we choose a frame where the incident spacelike photon carries q + = 0: q = [0. (qj_ + AjJ 2 /(P + + ((M 2 + 
A 2 L )/(l-C)P+,qJ andg' = [(P+ , (q ± + A ± ) 2 /(P+ ,q ± + A ± ]. 

In deeply virtual Compton scattering (DVCS) where Q 2 = —q 2 is large compared to the mass M and —t, one 
obtains Q 2 /(2P ■ q) = £, i.e. ( plays the role of the Bjorken variable in DVCS. For a fixed value of —t, the allowed 
range of £ is given by 

( -A / / A ]\/T2 \ 

(3) 

In the leading twist ignoring interactions at the quark-gauge boson(photon in this case) vertex, the amplitude 
contributing dominantly to Compton scattering in the deeply virtual region is given by 

M»" = M» v + M» v , (4) 

where M^ v and M£ u are the s(= (P + q) 2 )- and u{— (P — <7') 2 )-channel amplitudes (see Fig. [I]), respectively. The 
s-channel amplitude is given by 

M^ = -iN c e 2 q J -^Tr[ 75 (# + m) 7 "(#+ 4+m)>f(t- A + m) l5 {- P+ ff + m)] 

x H cov (k,P)K ov (k~A,P-A) 

[k 2 - m 2 + ie] [(k + q) 2 - m 2 + ie] [(k - A) 2 - m 2 + ie] [(P - k) 2 - m 2 + is]' 

where N c is the color factor and H cov (k, P)[H' cov (k — A, P — A)] is the covariant initial[final] state meson-quark 
vertex function that satisfies the BS equation. As usual in the LFBS formalism, we assume that the covariant vertex 
function H cov (k) does not alter the k~ pole structure in Eq. (jHJ) . The u-channel amplitude can be easily obtained by 
M£ v = AP s iU (q — » —q')- 

In deeply virtual q~ ~ Q 2 /(,P + limit, one obtains {jk+ jq + m) ~jq ~ 7 + <7~/2 in the trace term of Eq. (||). 
Consequently, combining this with the 2nd term of the denominator in Eq. (|5|) leads to q~ /\(k + q) 2 — m 2 + ie] ~ 
l/[P + (x — C + i £ )]i where x = k + / P + . Similarly, one can obtain (— q')~ /[(k — q') 2 — m 2 + ie] ~ l/[P + (x — ie)] from 
the u-channel amplitude. 

Adding these two s- and w-channel amplitudes, we obtain the Compton scattering amplitude in DVCS limit as 
follows 




M" = efr'M^ = -^e 2 / ™ ( — ^— + ^ ) H cov (k, P)H' cm {k - A, P A) 



iN c 2 r d 4 k ( 1 1 
2P+ Cq J (2tt) 4 \x- ( + ie + x~ie ; 
x Tr[ 75 (ff + m) f J '(J6- A + m) 75 (- P+ f6 + m)] 
[A: 2 - to 2 + ie] [(& - A) 2 - m 2 + ie] [(P - fc) 2 - m 2 + ie] ' 

For circularly polarized(e + = 0) initial and final photons^/, J are f or |), we obtain from Eq. 



(6) 



(ei-6l J )7 + + i(dxeI J ) 3 7 + Ti5, ( 7 ) 



where we use the identities ( 7 + ) 2 = 0, 7 1 7 + 7 1 = >y 2 <y+>y 2 = <y+ and 7 1 7 + 7 2 = — ^j+j 1 — i 7 + 7 5. Equation (0) 
reduces to 7 + (l ±75) for the parallel helicities(i.e. + for and — for J, J.) and zero otherwise. Since the axial current 



2 As discussed in |g], for a longitudinally polarized initial photon, the Compton amplitude is of order 1/Q and thus vanishes 
in the limit Q 2 — * 00. 
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7 + 75 does not contribute to the integral, i.e. ~ (kj_ x AjJ after the trace calculation, we shall omit this term from 
now on. 

Then, the DVCS amplitude(i.e. photon helicity amplitude) can be rewritten as the factorized form of hard and soft 
amplitude 

M^(P,q,P') = M^(P,q,P') = -e 2 q [ dx( ±— + —L-) ^(C, x, t), (8) 

H J \ x — Q + i£ x — is) 

where 

72 



^(C,M) = ^^H cov (k,P)H' cov (k-A,P-A) 

Tr[ 75 (^ + m)7+(^- ^ + m)j 5 {- P+ g + m)] 
[fc 2 - m 2 + ie][(k - A) 2 - m 2 + ie][(P - k) 2 - m 2 + ie] 



(9) 



The function J r „^,x,t) is so called "generalized parton distributions" and it manifests characteristics of the ordi- 
nary(forward) quark distribution in the limit of ( — > and t — > 0. On the other hand, the first moment of the 
•FjrCCi x : i s related to the form factor by the following sum rules |,§: 

dx 

T n (C,x,t)=F 7I (t). (10) 



o 1-f 

In general, the polynomiality conditions for the moments of the GPDs |^3|,Q defined by 

dx 



1 - 



x n - 1 T 7I ((,x,t)=F n ((,t) (11) 



require that the highest power of ( in the polynomial expression of F n {C,,t) should not be larger than n. These 
polynomiality conditions are fundamental properties of the GPDs which follow from the Lorentz invariance. We 
comment on how our model calculations satisfy the polynomiality conditions in Section IV(numcrical results). 

An important feature of the DVCS amplitude given by Eq. (||) is that it depends only on the skewedness parameter 
C = Q 2 / (2P-q) for large Q 2 and fixed it |(< Q 2 ), i.e. DVCS is equivalent to an inclusive process exhibiting the Bjorken 
scaling Xgj = C Note also from Eq. (H) that the imaginary part of the DVCS amplitude is proportional to C t). 

The single spin asymmetry(SSA) fl!Jg(| that can be measured in the scattering of a longitudinally polarized probe on 
an unpolarized target is proportional to the imaginary part of the amplitude, i.e. the value of JyfC; C t)- The recent 
measurements of SSA for the proton target have been reported by HERMES and CLAS J26| collaborations. We 
discuss the value of !F W (C, 1) m Section III C. 



III. GENERALIZED PARTON DISTRIBUTIONS IN LIGHT-FRONT BETHE-SALPETER APPROACH 

The GPDs entering as nonperturbative information in DVCS as shown in Fig. |l| can be represented by Fig. ||(a) 
where the small white blob shown in Fig. ||(a) represents the composite(nonlocal) operator || at the quark-gauge 
boson vertex. Since the longitudinal component, A + = QP + , of the momentum transfer A = P — P' is in general 
nonzero, the covariant diagram Fig. ^|(a) is represented by the sum of the LF valence diagram (b) defined in £ < x < 1 
(DGLAP) region and the nonvalence diagram (c) defined in < x < £ (ERBL) region. As one can see from Fig. ||(b) 
and (c), the large white blobs at the meson-quark vertices in (b) and (c) represent the ordinary LF wave function. 
However, the large black blob in (c) cannot be represented by the ordinary LF wave function since it is no longer a 
bound state vertex but non-wave-function vertex. This non-wave-function vertex causes the main source of difficulty 
in representing the GPDs(as well as the timelike form factor) in terms of light-front wave function. However, in a 
covariant BS formalism it can be represented as an analytic continuation of the usual BS amplitude. 



Note that our definition of ^-(Ci x i t) in this work is different from the one used in our previous work jjj by a factor (1 — C/2) 
in the normalization. We also used the notation £ for the "skewedness" parameter instead of previous £ to be consistent with 
the notation used by Radyushkin ||. 
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FIG. 2. Diagrams for GPDs in different kinematic regions for the case £ > 0: The covariant diagram (a) corresponds 
to the sum of the LF valence diagram (b) defined in DGLAP(£ < x < 1) region and the nonvalence diagram (c) defined in 
ERBL(0 < x < C) region. The large white and black blobs at the meson-quark vertices in (b) and (c) represent the ordinary LF 
wave function and the nonvalence wave function vertices, respectively. The small white blob at the quark-gauge boson vertex 
indicates the nonlocality of the vertex. 

In our previous work [12| , we have derived the GPDs of the pion starting from the covariant BS amplitude of the 
current ^(see Eq. (12) in which is essentially the same calculation of a;, i) given by Eq. ([)]). Since the 

detailed procedures for obtaining the effective solution for the non- wave- function vertex have been given in fl2]| , here 
we only briefly present the salient points of our previous effective method |P2] ] before we discuss the new treatment 
of the quark-gauge boson vertex function to ensure the continuity of the GPDs at the crossover(x = £) between the 
DGLAP and ERBL regions. 



A. Brief review of our previous method |12] 



The essential feature of our approach is to consider the light-front wave function as the solution of light-front 
Bethe-Salpeter equation(LFBSE) given by 0,f6|J|7|j2|] 



(M 2 - A^V(^,k' J± ) = / 1 [rfy'][d 2 l' ± ]/C(^,k' J± ;y;,lV)x'(y;,l' J ± 

Jo 



(12) 



k' 2 L )/x'- 



where^j JC is the BS kernel which in principle includes all the higher Fock-state contributions, Mq 2 = (to t jv j_ i 
(to 2 + 'k! 2 L )/{\ — x') is the invariant mass, and x'( x 'n k'i_i_) is the (final state) BS amplitude with the internal momenta 
of the (struck) quark for the final state, x' = (x — C)/(l — C) an d = k^ + x'Aj_- The internal momenta of the 
(struck) quark after the kernel are given by y 1 = (y — 0/(1 — C) an d = + so that the integration of y' runs 
from to l(or y runs from £ to 1). Note also d 2 Y ± = d 2 l± for a given y' . We define the valence BS amplitude(i.e. 
x > C) as X(2^2) an d the nonvalence BS amplitude(i.e. x < C) as X(i^3) where the subscript indicates the parton 
number before and after the kernel. Both the valence and nonvalence BS amplitudes are solutions to Eq. (|1J). As 
illustrated in Fig. ^(c), the nonvalence BS amplitude is an analytic continuation of the valence BS amplitude. In the 
LFQM the relationship between the BS amplitudes in the two regions is given by [p^ : |l(i| 



(M 2 -M^ 2 )x' (1 ^ 3) (^,k' l± ) = / [dy'][d 2 U]/C(^,k' i± ;^,lV)X(2^2)(y^l' J ±) J 



(13) 



where again the kernel includes in principle all the higher Fock-state contributions because all the higher Fock 
components of the bound-state are ultimately related to the lowest Fock component with the use of kernel as illustrated 
in Fig. |. 



4 We present the LFBSE for the final state meson-quark vertex with the final state momentum variables. However, the same 
form of LFBSE can be used in the initial state meson-quark vertex. 
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l < K 

v (l->3) /l (2->2) 

FIG. 3. Non- wave-function vertex(black blob) linked to an ordinary light- front wave function (white blob). 

After the Cauchy integration over k~ in Eq. (||), we obtain the valence(£ < x < 1) and nonvalence(0 < x < £) 
contributions to the GPDs of the pion as follows 



and 



^r(CM) 



Nr. 



2x(l — x)x 



/^i 2 k 
^JX(2^2) [X, k±)5+ i (X, k_L)x' (2 ^2) {?', k l), 



(14) 



-Nr. 



2x(l-x)(l-x") 
1 dy> 



16tt 3 



o 2/(1-2/) 



d 2 i x )C(x',-k!r,y',i , x )x' (2 ^ 2) (y',i'±), 



(15) 



where Eq. (13) has been used for the non- wave-function vertex as shown in Fig. [|[or Fig. ^[(c)]. The trace terms of 
the quark propagators, i.e. S^j(x,kj_) and S^ v (x, kj_), are given by Eqs. (21) and (25) in Ref. ]12j| , respectively. The 
light-front vertex functions of the initial [final] hadron X(2— >2)[X(2— >2)] an0 - the S au S e boson x 9 are given by^] 



X(2->2)(»jk±) = 



l LF 



M 2 



M 2 



LF 



M 2 



X ff (x",k" ± ) = 



(16) 



where x" = x/(,k." ± = k±+x"A±, and Mq 2 = (k"^_ + m 2 )/x" + (k"^ + 77i 2 )/(l — x") is the invariant mass of qq pair 
annihilated into the gauge boson. Equations (jlj) and ( |l5| ) were our essential results from our previous analysis [|12) . 
bmce the Eqs. @ and @ are divergent themselves with a constant LF vertex function /ilf, our idea is to replace 
the /iL^[or equivalently X(2^2) ( x , kj_)] with the standard LF vertex function(i.e. dressed vertex) |{l^[l4]|lqjl^]29fl , 
which has been successful in predicting many static properties of ground state mesons. However, we left the gauge 
boson vertex function x 9 as the bare vertex given by Eq.(|l6|) in our previous work [[L^ll . We(as well as the authors 
in PH) attributed this different treatment between gauge boson and meson vertex functions to the reason for the 
artifact of the discontinuity in the GPDs at the crossover(x = £) between the DGLAP and ERBL regions. 



B. Consistent treatment of gauge boson and meson vertex functions in LFQM: Update 



As discussed in our previous work |L2l and the works by others [ pO|JlC|Jl 1[| , in principle one should consider the same 
BS kernel K. for qq pair annihilating into external gauge boson as that for qq pair merging into the final state meson 
given by Eq. dl5|). This is illustrated in Fig. [| and we shall denote the GPDs coming from Figs. f|(a) and (b) as 
■^""(a) (£' x > (gi ven by Eq- (|l5|)) and .T 7 ™^ (£, X, t), respectively. Adding both contributions of T™^ and T^? b \ assures 
the cancellation of any infrared divergence that might occur in the kernel /C. 



5 Note that we slightly modified the definition of \ 9 m present work from the one in [[L2|, i.e. we take the common factor 
C/(l — C) from Eq. (24) in Ref. jL2| and put back into the prefactor leading to the present Eq. (|l5|). 



G 



i i-c i 

(a) (b) 
FIG. 4. Possible light-front time-ordered nonvalence diagrams with the kernel for qq pair merging into the final state meson(a) 

and annihilating into external gauge boson(b). 

The GPDs from Figs. |(a) and (b) can be written generically(dropping, for simplicity, relative label for internal 
momenta) 

£dxT^ a) (( : x : t) =J Q C dx J d 2 k ± j\y J d 2 l ±X f {v, U)M*, k±5 V, U)x ff (x, k ± )S+ v (x, k ± )^ a (x, k ± ), (17) 



dy^ b) (C,y,t) = J q dy J dH x j\x J d 2 k ±x { (x, k x )K b (y, U; x, k ± )x 9 (y, 1±)S+ (x, k x )x\(x, k ± ), 



(18) 



where the superscript i(f) in x indicates the initial(fmal) valence light-front wave function(i.e. white blob in Fig. [|) 
and the notation of x\(h) implies the product of intitial x\n,) m -^^S- @ a (b) and other terms(such as prefactor in 
Eq. (|l5|)) that depend on the internal momentum variables. The kernel ICb{y, lj_; x, kj_) in Eq. (|l8| ) is equivalent 
to /C a (x, kj_;2/, in Eq. ([Tt]) upon exchange of dummy variables (x,kj_) <-> (y, except an overall sign due 
to an exchange of quark and antiquark. If, for instance, the kernel is approximated by one boson exchange, i.e. 
K-a = g 2 9(y-x)/(l+ -k+)[(P-k)- -(l-k)~ -(P-l)~] and K b = -g 2 9(x - y)/(k + — l + )[k~ —l~ — (k — l)~], where 
y = l + 1 ' P + and x = k + / P + , then one can easily see K. a {x,y) = —fCb{y — ► x,x — ► y). Note also in this simple one 
boson exchange case that there is only one light-front time ordered diagram in each kernel and other contributions 
such as 9(x — y)-term in IC a and 6{y — x)-term in K, b vanish since they contribute to pair creations from the vacuum. 
Likewise, the vertex functions x f and x 9 m Eq. (|l^ ) are the same as those in Eq. ( [l8|) upon exchange of dummy 
variables (x,kj_) «-> (y, lj_). 



Exchanging (x, kj_) <-> (y, l±) in Eq. (|l_S 



d 2 k. 



dy 



), we now combine these two contributions, T n y > and T n y,s, to get 
d 2 l±X f (y, k ± ; y, 1_l)x 9 (x, kjJS+^x, k^x^x, k x) 



C 



1 - 



5^(x,k ± ) xLO^k-O 



(19) 



where /C = /C a and the relative(— ) sign between the two contributions assures the removal of infrared singularity 
at x = y and k^ = Ij_ (sec also Ref. [(L0|). In Eq. ([l9|), the (x,kj_) dependence in T™^ comes solely from the 
term K,{x, k±; y, l±)x 9 {x, kj_) (i.e. upper small loop in Fig. |](b)), which in principle includes the sum of all possible 
intermediate mesons being coupled to the external current. 

By putting the relative label of internal momenta back into Eq. ( pj| ) , we obtain the same form of J 7 ™"" as that given 
by Eq. © 



*r(c,z,*) 



-N c 



2x(l 
l 



-x)(l 
dy' 



d 2 k. 



x") J 16tt 3 



X(2^2)(x, k ± )S+ (x, k ± ) X g (x", k" ± ) 



y'{i-y') 



d 2 \ x £{x',k!r,y',W)x'(2^){y'Jx), 



but now with 



lC(x',k' ± ;y',l'±) =/C(x',k' ± ;y',l'j 



S£ v (v,i±)x{(y,i±) 



S' r L(x,k_ L )xj l (x ! kj 



(20) 



(21) 



7 



Thus, in order to make our previous model |12| more complete, the kernel K,(x', k'j_; y', V±) given by Eq. (||) has to 
be replaced by IC(x' , k'j_; y' , which is free from any infrared singularity that might occur in the kernel JC. Since it 
would be a formidable task to solve the kernel K. directly, we follow the technique illustrated in our previous analysis 
|l2| approximating G n defined by 

Grr = f ,,f y ' / d 2 l ± JC(x\k' ± ;y',l'±)X(2^)(y',l'±) (22) 

Jo yV--y)J 

as a constant. As we did in our previous analysis |T^], we discuss the validity of a constant approximation in the 
next section (Section IV) where our numerical analysis is presented. Although the technical aspect of our numerical 
analysis remains same, the replacement of kernel from AC to AC amounts to the addition of new contribution and 
thus effectively allow us to change the gauge boson wave function x 9 (a;",k"j_) from the simple energy denominator 
given by Eq. ([l6]) to the LF vertex function identical to the meson vertex function without introducing any additional 
parameters [| 

Comparing X{2— >2) with our light-front wave function given by Ref. |i"3y , we identify 



, , , /87T 3 fdkz [X(1-X)]V 2 lf , , 

where the Jacobian of the variable tranformation k = (k z , kj_) — > (x, kj_) is obtained as dk z /dx = Mq/[4x(1 — x)] and 
the radial wave function is given by <fi(k 2 ) = ^/ 1 /7r 3 / 2 /3 3 exp(— k 2 /2/3 2 ), which is normalized as J d 3 k\<fi(k 2 )\ 2 = 1. 
Note that the radial wave function 0(k 2 ) is essentially the same as the Brodsky- Huang-Lepage [[LLj LF wave function 
4> BUL = exp(— Mq/8/3 2 ) up to a constant factor. 

Thus, the effective gauge boson wave function x 9 ( x " 1 k"j_) is given by 



N r _ V dx" Mq 



where dk"/dx" = Mq'/[4x"(1 — x")\. Here, the radial wave function tfi 9 (x", k"j_) is given by 



'CO = \l ^33 ex P (-k" 2 /2/3 2 ) exp(A 2 /8/3 2 ), (25) 



where k" 2 = k" 2 + k'fi and the virtuality A 2 = t ^ 



Substituting Eqs. © and @ into Eqs. Q and @, we obtain the valence and nonvalence contributions to the 
GPDs of the pion in LFQM as follows 



y k]_ + mVk^ + m- 2 




i,kiy(i",k*i) 

x ^-K^^ + xjl x)x> { M 2 -M 2 ) f dyl r H [KS^MM mV±)i (27) 



V'k^ +m 2 v / k" 2 L + to 2 yT'l + to 2 



where Eq. (|2^) is the same as our previous formula Eq. (28) in |Q but Eq. ( f27j ) is now different from Eq. (29) in |§ 
due to the modification of the gauge boson vertex function. It is worthwhile to note that the nonvalence T™ in 



6 In principle, one should use the vector(such as p) meson wave function replacing the gauge boson vertex function. We do 
not need to introduce any additional parameters since we used the same gaussian wave function with the same parameter j3 
for both 7r and p mesons in our previous LFQM analysis |13fl. 
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Eq. ( p7| ) receives not only the on-mass shell propagating contribution (i.e. the term proportional to (kj_ • k'j_ + m 2 )) 
but also the instantaneous contribution(i.e. the term proportional to x(l — x)x'(M 2 — Mq )) from the spectator quark 
propagator. This also contrasts to the valence T™ 1 given by Eq. (p6|), which receives only on-mass shell propagating 
contribution. As we shall show in our numerical calculation, the instantaneous contributions from the nonvalence T™ 
becomes substantial for large C > 0.5. 

As discussed earlier, we shall treat the last dy'd 2 lj^-integral term in Eq. ( f27|) as a constant G^ and check if indeed 
a constant G-n approximation is valid by varying the value of £ in the sum rule expressed in terms of J-% al and T™ , 
i.e. 

FAt)= f -^ T Mi,x,t)= f JE T jy*(£ )X)t )+ f ( -^fr((,x,t), (28) 
Jo 1 — 5 J( 1 — I Jo 1 — | 

for given — t. Varying the value of £ in this sum rule, Eq. (^8|), we can also check the frame-independence of our model 
as we discuss in the next section(Section IV). 

C. The GPD value at the crossover: 7^(^,(^,4) 

Although the continuity of GPDs at x — £ is assured in our formulation presented in the last subsection (Section 
III B), the value of ^(C, C> t) vanishes in our model calculation as we shall see in our numerical results (Section IV). 
The reason why this occurs in our model calculation for the meson target is because the final state meson (depicted 
in Fig.^(b) and Fig.||) shares the total longitudinal momentum fraction 1 — C only between the two constituents so 
that, as x — £ — > (i.e. the struck constituent loses its longitudinal momentum completely), the single spectator- 
constituent carries the total longitudinal momentum fraction 1 — C an d the final state two-body wavefunction at this 
kinematical point (x — C) is zero. However, the situation is entirely different for the three-body wavefunction such 
as in the proton target because even at x = Q the two spectator-constituents can still share the total longitudinal 
momentum fraction 1 — £ between themselves and the final state three-body wavefunction does not vanish at x = £. 
As discussed in the literature [|l9|,^0) and mentioned in Section II, the value of GPDs at x = C is directly related to the 
SSA that can be measured in the scattering on an unpolarized target by the lepton polarized parallel or antiparallcl 
to its direction. The recent measurements by HERMES and CLAS collaborations showed that the SSA is indeed 
non-zero for the proton target and our observations based on the three-body wavefunction are not inconsistent with 
these experimental evidences for the non-zero GPDs of proton at x — £. 

By the same reasoning, if one includes the quark-antiquark-glue three-body state (and any other multi-constituent- 
states) beyond the ordinary two-body (QQ) constituents in our model calculation, then it is in principle possible to 
get a non-zero .^(C, £, i). I n the chiral-quark-soliton model it is shown that JvCC Cj^) is non-zero. Thus, if 
the experiment on the pion is possible at all, then the SSA measurement would be crucial to discreminate different 
models and also find how much percentage of multi-constituent-component is necessary beyond our current QQ state 
for the more realistic model construction. However, it is interesting to note that our current model calculation still 
satisfies the polynomiality conditions (See Eq.(|ll])) and our result on the isosinglet GPD of the pion is qualitatively 
very similar to the result obtained by the chiral-quark-soliton model satisfying the soft pion theorem, as we comment 
in the following Section (Section IV). 

IV. NUMERICAL RESULTS 

In our numerical calculations, we use the model parameters (m,j3) = (0.22,0.3659) [GeV] obtained in Ref. |l3| for 
the linear confining potential model. 

In Fig. we show the comparison of the ^-dependence of G T between the two gauge boson wave functions, i.e. 
our new gaussian type (black data), 4> 9 ~ exp[(i — Mg 2 )/8/3 2 ], and the previous [Of monopole type (white data), 
ipa l/(t — Mq 2 ) (see Eq.(|l6|)), for a few different momentum transfers, — t = (circle), 0.5 (square), and 1.0 (triangle) 
[GeV 2 ], respectively. The invariant mass of the go pair at the gauge boson vertex is represented as M 2 = Mg' 2 (see 
Eq. (|16|)) in the figure. As one can see in Fig. g, both (black and white) G^ values show approximately constant 
behavior for small momentum transfer —t<l GeV 2 region. It is not surprising to see that G w becomes very large as 
C — > 0, however, this does not cause a significant error in our calculation because the nonvalence contribution in the 
very small £ region is highly suppressed. Also, a rather large fluctuation of black data for a rather large momentum 
transfer region — t = 1 GeV 2 (even exhibiting a sign change near £ = 0.3) is neither unexpected nor troublesome 
because the nonvalence contribution is significantly suppressed in this large momentum transfer region where the 
gaussian wave function is also exponentially reduced. In order to check the reliability (i.e. frame-independence) of 
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our constant approximation, we compare the exact solution of the pion electromagnetic form factor obtained from 
£ = value with that of nonzero £ value using our constant (average) G n approximation. 



yj 1 

CD 



• For fixed -t=0.0 GeV 
■ For fixed -t=0.5 GeV 2 
a For fixed -t=1.0 GeV 2 



O O O 



6 a o 
o * 



~1/(t-M 



n m A B . W B a 

■: iii**** 



~exp[(t-M 2 J8 p 2 ] 
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0.2 



0.4 



0.6 



1.0 
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5. Comparison of the ("-dependence of between our new gaussian type(black data) and previous monopolc 
(white data) gauge boson wave functions for different momentum transfers — t — (circle), 0.5(square), and 1 (triangle) 



[GeV ], respectively. The invariant mass of the qq pair at the gauge boson vertex is represented as M, 



go 



Mr! 



Q, V! f0.12 



0.2 



0.0 



0.0 



Exact Sol. K=0) 

Effective Sol.(full) with (=0.3 

Effective Sol.(full) with (=0.6 

o Amendoliaetal. (1986) 
• JLab (2001) 
Nonvalence with (=0.6 




0.2 



0.4 0.6 
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FIG. 6. The effective solution of pion form factor using a common average — Grave. = 0.12 value with gaussian type gauge 
boson wave function for ( = 0.3 (thin solid line) 0.6 (long-dashed line) compared with the exact solution (thick solid line) as 
well as experimental data p2j,B3j. The dotted line represents the nonvalence contributions to the form factor for ( = 0.6 case. 



In Fig. O, we show our effective solutions of the pion form factor with gaussian (f> 9 (x", k"j_) for £ = 0.3 (thin solid 
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line) and 0.6 (long-dashed line) cases obtained from our average value of G^ = G avo . = 0.12 and compare them with 
the exact solution (thick solid line) as well as the experimental data j3^j3j| . The dotted line represents the nonvalence 
contributions to the form factor for £ = 0.6 case, which increase (decrease) as £ gets larger (smaller). Also, there 
are — t m i„ values for nonzero £ due to A^ > (see Eq. @). We thus use the analytic continuation by changing A_l 
to iA± in Eqs. ( ptf ) and (27) to obtain the result for < — t < — t min where there is no singularity. A continuous 
behavior of the form factor near — i m j n confirms the analyticity of our model calculation. As far as the G w and the 
form factor calculations are concerned, the results obtained from the gaussian type gauge boson wave function are 
overall not much different from our previous results fl2]| obtained from the monopole type gauge boson wave function. 
However, as we shall see below, they are distinguished by the calculation of the GPDs. 




FIG. 7. Generalized parton distributions of the pion at — t = 0.5 GeV 2 with £ = 0.3 in (a), and 0.6 in (b), respectively. 

In Fig. |7[ we compare the nonvalence contributions (in ERBL region) to the GPDs J-^{C, x, i)/(l — C/2) of the pion 
obtained from the gaussian type gauge boson wave function(thin solid line) with those obtained from the monopole 
type wave function(long-dashed line) for fixed momentum transfer —t = 0.5 GeV 2 but with different skewedness 
parameters ( — 0.3 in (a) and 0.6 in (b), respectively. The thick solid line represents the valence contribution(in 
DGLAP region), which is common for both gaussian and monopole type gauge boson wave functions. Our results are 
obtained from average values of G T , i.e. G^ = 0.12 for the gaussian type gauge boson wave function and 0.32 Jl| for 
the monopole type gauge boson wave function. We also plot the instantaneous contributions to the GPDs with the 
gaussian type(dotted line) and the monopole type(dot-dashed line) gauge boson wave functions for ( = 0.6 case. We 
note that the results of the GPDs shown in Fig. |?] with the average G^ values for both gaussian and monopole type 
gauge boson wave function cases are very close to the exact ones, i.e. frame-independent, as one can deduce from 
Fig. || for the gaussian wave function and Fig. 9 in []l2| for the monopole type wave function, respectively. 

As shown in Fig. |?], the most dramatic change in our updated result of the pion GPDs is that the solutions with 
the gaussian type gauge boson wave function show the continuity at the crossover(x = £), while the solutions with 
the monopole type |l2[ gauge boson wave function show discontinuity at the crossover. As we discussed before, such 
discontinuity at x = £ with the monopole type gauge boson wave function is just an artifact due to the different 
x — > C behavior between the gauge boson vertex (x 9 ) and the hadronic vertex (x(2->2)) functions (see Eq.(|l6|)). This 
discontinuity problem of the GPDs appeared in our previous analysis [fl2"[ is now resolved by realizing that the same 
type of bound state wave function should be used for both hadron and gauge boson. Furthermore, we can obtain 
finite results for the scattering amplitudes for the DVCS region given by Eq. (||) since both valence and nonvalence 
solutions of the GPDs are not only continuous but also vanish at the crossover^ = £) as well as at x = 0. As we 
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discussed in Section III C, ^(C, Cj t) can be in principle non-zero if the multi-constituent-states are included beyond 
the ordinary QQ state. Within our current QQ model, however, we have checked the nth moment F n (£, t) in Eq.(|ll|) 
for up to n = 3 as well as the isosinglet GPDs. With our updated formulation presented in Section III B, we did not 
find any difference from our previous results, Figs. 10 and 11 of Ref. []l2| , based on the formulation summarized in 
Section III A except the change of value discussed above. Thus, we do not show those figures in duplication but 
have confirmed that the polynomiality conditions are satisfied in our current model and our result on the isosinglet 
GPD of the pion is qualitatively very similar to the result (Fig. 5 of Ref. plj ) obtained by the chiral-quark-soliton 
model including the D-term in the ERBL region. The D-term generates the highest power in the polynomial, e.g. 
the second moment of the pion isosinglet GPD is given by F2 =0 ((,t = 0) = |(1 + C( 2 ), where C — —1/4 in the 
chiral limit 00 Our model calculation from Figs. 10 and 11 of Ref. Jl| gives C = —0.2843. Incidentally, we note 
from Fig. 5 of Ref. [ pT| that the chiral-quark-soliton model not neglecting the dynamical quark mass generated in the 
spontaneous breaking of chiral symmetry also gives the value of |C| slightly greater than 1/4. 



3. i 
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(a) 



IM lot l with <t> ~exp[(t-M j 



0.4 0.6 
-t[GeV 2 ] 




0.4 0.6 
-t[GeV 2 ] 



FIG. 8. (a) The real part of the DVCS amplitude M IJ = M TT + M u (see Eq. (||)) obtained from the gaussian type gauge 
boson wave function and the average = 0.12 for fixed C, = 0.6: The thick solid, dashed, dot-dashed lines represent the 
total squared amplitude Re(M 7J ) 2 , valence part of M IJ = M™i, an d nonvalence part of M IJ — M^v, respectively, (b) The £ 
dependence of the DVCS squared amplitude Ke(\M IJ \ 2 ). 

In Fig. |(a), we show the real part of the DVCS squared amplitude (thick solid line), Re(|M /J | 2 )=Re(|M TT +M u | 2 ) 
=4Re(|MTT| 2 ) in Eq. (§), obtained from the gaussian type gauge boson wave function and the average = 0.12 for 
the range of < — t < 1 [GeV 2 ] with fixed £ = 0.6. The thick dashed and dot-dashed lines represent the valence 
{M IJ = M' J j) and nonvalence (M IJ = M^) contributions to the total amplitude, respectively. It is interesting 
to note that the valence and nonvalence amplitudes interfere with each other destructively not constructively. We 
also note that the DVCS squared amplitude Re(|M /J | 2 ) for the monopole type gauge boson wave function shows the 
logarithmic divergent behavior as expected due to the discontinuity at the crossover x = £. We thus do not present 
the corresponding result for the monopole type gauge boson wave function. In Fig. ||(b), we show the £ dependence of 
Re(|M /J | 2 ) for the gaussian type gauge boson wave function. In our model calculations, the DVCS squared amplitude 



7 The skewedness parameter £ in Ref. pi[ is related to £ by £ 



Note, however, that the difference between £ and £ is 



negligible in the limit t — > because the physical region of £ is restricted by Eq.(E 
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increases as £ decreases. As in the case of Fig. [| we confirm the analyticity of our calculation from the continuous 
behavior of our results near —t = — tmin- 



V. CONCLUSION 



In this work, we reinvestigated our previous light-front quark model analysis |i2fl of the GPDs in the deeply virtual 
Compoton scattering (DVCS) of the pion, 7*7r — > 771- . We improved our previous effective treatment [Q of handling 
nonvalence (or higher Fock-states) contributions to the GPDs with the inclusion of the BS kernel for the qq pair 
annihilating into the gauge boson(see also _pO| , ^0 11 for a similar discussion) in addition to that for the qq pair 
merging into the final state meson(see Fig.^^It leads us to have a vector meson dominance picture at the gauge 
boson vertex within our constant approximation of G v in Eq. ( p2| ) and in turn to use the same type of gaussian 
wave function at both gauge boson and meson vertices. The continuity of the GPDs at the crossover(x = Q between 
the DGLAP and ERBL regions is ensured by a consistent treatment of the gauge boson and meson wave functions. 
Subsequently, the real part of the DVCS amplitude M IJ is calculated and its finiteness is confirmed (see FigsJ|(a) 
and (b)). In our current QQ constituent model, we obtain ^-(CjCi^) = indicating that the single spectator cannot 
share its longitudinal momentum at x = Q. In principle, J-^iCi C? t) can be non-zero if our model is extended to include 
the multi-constituent-components. However, our current model is capable of satisfying the polynomiality conditions 
and yields qualitatively very similar results on the isosinglet GPD of the pion indicating that the D-term is effectively 
included in the ERBL region. Also, one should note that our model is not inconsistent with the recent HERMES and 
CLAS data on the SSA for the proton target. 

Finally, the GPD of the pion is invariant under the time reversal symmetry(A — -> — A)^|. As illustrated schematically 
in Fig. |9|, it is not difficult to see that our effective method does not violate the time reversal symmetry. 




FIG. 9. Illustration of time reversal(A — > —A) process of the GPDs in Fig 




Using our model calculation of the DVCS amplitude in which the final photon is emitted by the pion, one can 
calculate the differential cross section of the virtual Compton scattering of the pion, i.e. ne —> 7TC7, by incorporating 
the Bethe-Heitler(BH) process f34j,||j| in which the final photon is emitted by either the incoming electron or the 
outgoing electron. Considerations along this line are in progress. 
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